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Abstract 

We examine the bubble nucleation rate in a first-order phase transition 
taking place in a background Jordan-Brans-Dicke cosmology, We compute 
the leading order terms in the nucleation rate when the Jordan-Brans-Dicke 
field is large (i.e., late times)^by means of a Weyl rescaling of the fields in 
the theory. We find, that despite the fact that the Jordan-Brans-Dicke field 
(hence the effective gravitational constant) has a time dependence in the false 

vacuum, at late times the nucleation rate is time independent. 
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Interest in theories with variable gravitational constant has been rekindled recently 
with the advent of Kaluza-Klein [l] and string theories [2], Theories from extra dimen- 
sions typically have a scalar field, usually called the dilaton, coupled to the Ricci scalar 
in a way reminiscent of Jordan-Brans-Dicke (JBD) theories [3]. The JBD theories have 
often been used as toy models to gain insight into more complicated theories. One of 
the most obvious questions to ask is how these theories differ from standard Einstein 
gravity in their cosmological evolution. We expect that the cosmological behaviour of 
such theories at early times will in general be quite different from the standard cosmology 


[4]- 

Our main interest here involves the nucleation of bubbles of true vacuum in first-order 
phase transitions. Steinhardt and La [5] have recently used a JBD theory coupled to a 
scalar field theory admitting a false vacuum to try to save the old inflation model of Guth 
[6], Their model, dubbed extended inflation, represents an interesting new addition to 
the existing panoply of inflationary models [4] (despite the fact that it too may have some 
problems [7,8]), since the fine-tuning problems usually found in new inflation models [4,9] 
can be eliminated. 

In this letter we examine the effects upon the bubble nucleation rate of a time- 
dependent background JBD field. Although we only consider JBD theories, our results 
can be trivially extended to other models with varying gravitational constants, e.g., 
induced gravity models where the Ginzburg-Landau scalar field plays the role of the 
JBD field [10], or Kaluza-Klein and string theory cosmologies where the dilaton field 
plays the analogous role [11]. While one would expect that a complete calculation of the 
effects of the JBD field and the metric on the nucleation rate would yield a time dependent 
nucleation rate (due to the time dependence of the JBD field in the false vacuum), we find 
that at late times (or equivalently, for large values of the JBD field), the rate becomes 
approximately time independent. Thus, at least in this regime, the calculations of Refs. 
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[7,8,12,13] are applicable. 

Let us start with the action used by Steinhardt and La: 

/ _ Q <5>Q <J> 1 

c PxV=9 + + L g >“'d^d v v-V{ < T) . (1) 

Here <r is the scalar matter field that plays the role of the “inflaton” in extended inflation 
models. We assume that c has a potential V(<r) with a metastable minimum at a — <tfv 
with V(<r = <tfv) = py ^ 0, and a true ground state at a = <To with V(cr 0 ) = 0. $ is the 
JBD field (with dimension mass 2 ), and g ^ is the metric tensor. Any time dependence in 
the JBD field will be manifest as a time dependence of the coefficient of i?, which serves 
as an “effective” gravitational constant: $ = 1/167 tGeff- The theory is characterized 
by a dimensionless parameter u>; in the limit u> — ► oo the theory is identical to general 
relativity. Present experimental limits on lj require u> > 500 [14]. 

We will be interested in cosmological solutions to the theory where the metric tensor 
and the fields depend only upon time, and the line element takes the Robertson-Walker 
form ds 2 = dt 2 — a 2 (t)dl 2 , where a(t) is the cosmological scale factor and dl 2 is the spatial 
metric 

dr 2 

dl 2 = j— + r 2 d6 2 + r 2 sin 2 6d<f> 2 , (2) 


where k = 1 , 0, or —1 for a closed, spatially flat, or open Universe. 

With the assumptions of a Robertson-Walker metric, <r — constant, and $ = $(<), 
the field equations become 



where a dot denotes djdt and py = —py = V(<r = <rpv) is the (constant) vacuum energy. 
We will assume that the curvature term in Eq. (3) can be neglected and set k = 0. The 
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resulting equations have a solution in which the scale factor a(t) increases as a (large) 
power of time, and the JBD field also increases with time [15]: 


$(<) = $( 0)(1 + £ 0 2 
o(0 = a(0)(l + £tr +1/J > 


where 


B 



a 


L«*(o)J 

(3 + 2w)(5 + 6u>)/12. 


(5) 


(This is not the most general k = 0 solution, although a large class of solutions approach 
this form at large times [7].) 

So long as o' = constant, a(t) and $(t) will increase in time. If there is a metastable 
minimum at a = <r F v, eventually bubbles of true vacuum with <r « <r 0 will be nucleated. 

It is this nucleation rate which we wish to study. 

The field theoretic calculation of the decay of the false vacuum in flat space-time 
was essentially solved by Coleman and Callan and Coleman [16]. Incorporation of the 
gravitational effects upon the nucleation rate was studied by Coleman and DeLuccia [17]. 
However, the latter analysis is not directly applicable here because of the nonstandard 
gravitational action whose effective coupling constant changes with time. In order to 
sidestep these difficulties, let us recast the theory in a form with a time-independent 
gravitational coupling constant. To do this, we perform a Weyl rescaling [18] to remove 
the coupling of $ to R. In more detail, we define new fields g and V’ by 


9 nv — A (*) 9v» 

= *,ln[*(3 + 2u)/itf] , (6) 

where rp* = (3 + 2 u>)/ 16 tt G n with G N the presently observed value of Newton’s constant. 
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Choosing the conformal factor ft(i) to satisfy = l/167rC?^ changes the action to 

S[g,il),<T\ = J d*x R 

+ exp(—ip / ip Q ) l -r v d»<r d v cr - exp(-2 rfr/fo) V(a)j . (7) 

Thus, we now have a system consisting of standard Einstein gravity (with constant 
gravitational constant), a minimally coupled field and the inflaton field <r now coupled 
to if>. A similar approach was used by Gross and Perry [19] in their analysis of monopoles 
in Kaluza-Klein theories. 

How can a theory with a variable G become a theory with a constant G1 The simple 
answer is that there are two mass scales, a gravitational mass scale and a mass scale in 
but the only relevant observable is the ratio of the two mass scales. In the theory 
as written as in Eq. (1), the a mass is constant in time while the gravitational mass 
changes, while in the theory as expressed in Eq. (7), the gravitational mass is constant 
but the tr mass vanes. Of course, in both pictures the ratio of the gravitational mass to 
the <7 mass evolves in time in the same way, and one may choose to work with whichever 
picture is more convenient. 

Again making the assumption that the metric is of the Robertson- Walker form, <r = 
constant, and V> only depends upon time, the equations of motion are 

+ & = — 3 — [2^ +ex P(- 2 V'/V’oWJ 

r/," + 3 1^' = A exp (— )K. (8) 

In Eq. (8) we have denoted the time and scale factor in the rescaled system of Eq. (7) as 

< and d(t), respectively, and prime denotes d/dt. The k = 0 equations of motion again 
have a power law solution: 


a(t) = d(0)(l + Ct) 3/i+u/3 

= V>(0) + y>oln(l + Ct), 
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( 9 ) 


where here, 


C = 


1 


Pv 


1/2 


= 2(167rG^*(0)] _1/, B. 


( 10 ) 


a$(0) L 24 kGs* 

At this point we would like to see how the two solutions are related to one another. 
We first note that the conformal factor in our Weyl rescaling is given by 


n(0 = [16;rG^(0r 1/2 = ^(l + B0- 1 . 


( 11 ) 


Since the two line elements ds 2 = dt 2 — a 2 (t )dl 2 and ds 2 = dt 2 — a 2 (t)dl 2 are related by 
da/da = fl, we must have dl = dl and 


dt 

dt 



( 12 ) 


Integrating the di/dt equation and fixing the constant of integration so that the a{t)ja{t ) 
equation is satisfied, we obtain 


(1 + Ct) = (1 + Bt) 2 . 


( 13 ) 


Having completed these preliminaries, let us turn to the the bubble nucleation rate. 
First, we note that the conformal relation between the original and the rescaled theories 
fixes the relationship between the nucleation rates in the two theories. In the latter 
theory the nucleation probability per physical volume per time is 
dP{t) 

~ dV(i)dt ’ 

where P{t) is the probability of bubble nucleation and V(t ) is the physical volume 
measured in the rescaled system at time t. The equivalent rate in the original theory is 
then 






dP{t) 

dV(t)dt 



j t A(f) = sr<(!) A(< ), 


( 15 ) 
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where we have used Eqs. (11) and (13) and the fact that P{t) = P[t (£)]. 

In the JBD theory there are two effects which might cause the bubble nucleation 
rate to deviate from its flat-space value. The first is the gravitational coupling, while 
the second is the time dependence of the JBD field. In the rescaled formulation of the 
theory the gravitational interaction has the standard form, and so we might expect direct 
gravitational effects to be similar to those in a standard theory. In order to isolate any 
time dependence coming directly from the JBD field, let us consider the theory of Eq. (7) 
in the Gs — ♦ 0 limit. In this limit the first, gravitational, term in the action should 
certainly be omitted. What is perhaps less obvious is that the second term, containing 
the kinetic energy of the JBD field, should also be neglected. To make this a bit clearer, 
consider replacing i>(x) by rp(x) = V , ( a: )/V'o- The JBD kinetic energy term now has the 
form jV’o 9 ftu Because V’o is proportional to GJf 1 ^ 3 , this term scales in precisely 

the same way as the Ricci term as Gn — » 0. 

We are thus left with the last two terms of Eq. (7). Only the matter field a remains as 
a dynamical quantity. Without their time derivative terms, the metric and the JBD field 
now enter essentially as nondynamical, externally specified quantities. Since gravitational 
effects can only be small if the bubble size is much less than the curvature of space, we 
take the metric to be flat, and write the action of this truncated theory as 


S[a] = Jd*x ^b(d^) 7 + b 3 V(cr) , 


(16) 


where we have defined b = exp(— ■0/V'o)- We now wish to calculate the nucleation rate 
in this truncated theory, or, more precisely, to determine its dependence on the value 
of b. With the gravitational interaction no longer present, we can use the methods of 
Ref. [16], although with some modifications necessitated by the nonstandard form of the 
kinetic energy term. The first step is to find the “bounce” <r(®); i.e., the minimum action 
solution of the Euclidean equations of motion containing a true vacuum region where 
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<7 « <r 0 but with a approaching its false vacuum value at infinity. The bubble nucleation 
rate per unit volume is then given by A = >lexp( — Se), where Se is the Euclidean action 
of the bounce while the prefactor A is given by 


A = 


det' [S£(5-)] 


det ^(o-fv)] 


-1/2 


n 


2t t) 


1/2 


( 17 ) 


Here, the prime denotes that the functional determinant is to be evaluated in the subspace 
orthogonal to the four translational zero modes. The are normalization factors for 
these zero modes, defined so that the properly normalized modes are C , " 1 / 2 9 p o’, fi = 1...4, 
so that 


C M = J d^d^a) 2 , 


(18) 


with no sum over fi implied. (With an 0(4)-symmetric bounce, the are all equal and 
can be simply related to the bounce action; we will not need to use this fact.) 

By a simple rescaling of coordinates, one finds that the bounce solution has the form 
c(x) = a(y/bx), where a(y ) is the bounce solution when 6=1. By substituting this into 
the action and then rescaling the integration variable, we see that Se is the same as for 
the 6=1 case; all of the 6 dependence lies in the prefactor. Performing the functional 
variation of the Euclidean action gives 


j4det = 


det' [5g(*)] 

-1/2 

det' [— 6<? 2 + b 2 V"(a)] 

det [^(o-Fv)] 


det [-bd 2 + b 2 V"((r FV )] 


- 1/2 


(19) 


To determine the 6 dependence of the determinants, we note that if $*(x) is an eigen- 
function of the operator — d 2 + V"(a) with eigenvalue A, then ^(x/hx) will be an eigen- 
function of the operator — bd 2 + b 2 V"[a) with eigenvalue 6 2 A; the same relationship holds 
for the eigenvalues associated with the false vacuum configuration. Now if we note that 
the primed determinant has four fewer eigenvalues than the unprimed one, we see that 
Adet contains a factor of [(6 2 ) -4 ] -1 / 2 = 6 4 ; i.e., .Adet = 6 4 Adet> where Adet is the deter- 
minant ratio 6=1. Finally, the 6 dependence of the can be obtained by substituting 
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our bounce solution into Eq. (18) and rescaling the integration variable. The result is 
that = b~ 1 C^ i , where again the hat indicates the 6 = 1 value. The prefactor A thus 
has a net factor of b 3 , so that A (£ ) = 6 2 Ao, where Ao is the nucleation rate for a normal 
scalar field theory with potential V'(c). Since 6 is a function of the time-dependent JBD 
field, we thus have a time-dependent nucleation rate in the rescaled theory. However, 
this time dependence disappears if we go back to the original theory. Eqs. (6) and (11), 
together with the definition of ipo, imply that 6 = fl 2 (t). Using Eq. (15), we find that 
the nucleation rate of the original theory is A(£) = A 0 , and is thus time independent in 
the approximation to which we are working. Indeed, this is just the result one would 
have obtained by omitting both gravitational and JBD effects from the original action of 

Eq» (!)• 

It may be a bit unsettling that the nucleation rate can be either time dependent or 
time independent, depending on the framework in which one chooses to work, but this 
is inevitable once we exploit the possibility to make a time-dependent rescaling of space- 
time. The ambiguity is eliminated if one asks questions which are phrased entirely in 
terms of physical lengths. Thus, if we calculate the probability that a bubble will nucleate 
within a space-time volume equal to the fourth power of the <r-Compton wavelength, the 
result will be the same in either formalism and, in the limit we are working, will be 
independent of time. 

Similarly, we should require that both frameworks give the same answer to the ques- 
tion of whether and, if so, when the true vacuum percolates. Some indication of this can 
be found by considering the dimensionless parameter e = X/H*(t) where H(t) = a(t)/a(t) 
[7,8,12]. In a de Sitter universe described by general relativity, percolation occurs if this 
parameter is larger than a critical value, tcRi lying in the range 1.1 x 10 -fl < ccr & 0.24 
[12]; the same is expected to be true in the large-w limit of JBD theory. To compare our 
two theories, we calculate i(t ) = Using our formulae above, and making use of 
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the relation between t and £, we find that 


e 



2u> 4- 1 \ 4 
2cj -f- 3 / 


=(()■ 


( 20 ) 


so that the two theories agree in the limit of large u>. 

Let us comment on the validity of our approximations. We have worked in the limit 
Gff —* 0. Although we have not yet calculated the first order corrections to this limit, we 
can use the results of bubble nucleation calculations in standard gravity as a guide. In 
those calculations, the terms of order Gn are small provided that the mass scales in the 
matter part of the action are much less than the Planck mass, and that the bubble size is 
small compared to the curvature of space. By analogy, we expect our approximation to be 
reliable when the effective Planck mass induced by the JBD field is much greater than the 
mass scales associated with the <r field. Since $ increases with time, our approximation 
should be best in the large-£ limit. 

We should also comment on our assumption that ip could be treated as constant during 
the bubble nucleation process. This would imply that $ was also constant. However, 
we have explicitly used the fact that both a and $ were time dependent in order to 
derive our result for A. The resolution of these seemingly contradictory assumptions is 
that the (imaginary-time) bounce configuration used in computing the tunnelling action 
is distinct from the (real-time) backround metric and JBD field configuration governing 
the evolution of the Universe. It is this latter configuration that remains time dependent. 
A similar situation occurs in old inflation nucleation rate calculations, in which the 
metric is frozen out of the bounce solution, yet the real-time Universe is still expanding 
exponentially. 

To conclude, we see that at least at large times, the influence of the time dependence 
of the JBD field does not give rise to any time dependence of the nucleation rate. It 
would be of great interest to understand how to calculate the full consequences of the 
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JBD field on this rate. This is being done by us [20], but there are some questions of 
principle, such as how to specify boundary conditions for the JBD field in the bounce 
configuration. We also note that if one chooses to work in the rescaled system where 
gravity is normal, microphysics is abnormal, and one must take care in calculation of 
bubble nucleation. 

Gravitational effects on vacuum decay in JBD theories have also been considered by 
Accetta and Roman elli [21]. 
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